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Abstract. We study the QCD phase structure in the three-flavor Nambu-Jona-Lasinio model, 
incorporating the chiral-diquark interplay due to the axial anomaly. We demonstrate that for a 
certain range of model parameters, the low temperature critical point predicted by a Ginzburg- 
Landau analysis appears in the phase diagram. In addition, we show that the axial anomaly presents 
a new scenario for a possible BEC-BCS crossover in the color-flavor locked phase of QCD. 
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INTRODUCTION 

The phases of QCD at finite temperature T and quark chemical potential \x are being 
actively studied. In particular, at low T and high ji a color superconducting (CSC) 
phase [1, 2] characterized by a diquark condensate (qq) is expected to appear owing 
to the attractive interaction between quarks, provided either by one-gluon exchange or 
by instantons. On the other hand, when the system at finite /i is heated, a transition to a 
quark-gluon plasma (QGP) takes place at a (pseudo-)critical temperature. 

Our understanding of the phase structure based on the lattice simulations is still 
immature due to the severe sign problem at finite /i. Therefore the analyses so far have 
relied mainly on specific models of QCD, such as the Nambu-Jona-Lasinio model [3, 4], 
the Polyakov-Nambu-Jona-Lasinio (PNJL) model [5], and etc. Some of these model 
indicated the possible existence of a critical point located at high T [6]. 

The interesting possibility of a second critical point at rather low temperature in the 
color-flavor locked (CFL) phase was recently predicted on the basis of general Ginzburg- 
Landau (GL) analysis [7]. Moreover, this critical point has proven to make a quark- 
hadron continuity possible [8, 9, 10]. In two-flavor QCD, similar critical points have 
been found in the NJL model [11], although their origin is different from axial anomaly. 

This work is aimed at locating this critical point in the -phase diagram using 

the phenomenological NJL model [12]. Starting with the three-flavor NJL model incor- 
porating the axial anomaly induced chiral-diquark interplay, we study the location of 
the new critical point and its dependence on the strength of the anomaly. We also ob- 
serve that the axial anomaly triggers a crossover between a Bose-Einstein condensed 
state (BEC) of diquark pairing and Bardeen-Cooper-Schrieffer (BCS) diquark pairing 
[13, 14, 15, 16, 17] in the CFL phase. 



INCORPORATING THE AXIAL ANOMALY IN NJL MODEL 



The Lagrangian of the NJL model with three-flavors consists of three terms: 

if = q(iy^ -m q + Mo)q + i^ (4) + ^ (6) , (D 

where q = (u, d, s) is the flavor triplet quark field, m q is a flavor symmetric quark mass 

(m u = 111(1 = m s ), and ji is the chemical potential for conserved quark number, if ( 4 ) and 

^(6) 

are the four-fermion and six-fermion interactions, respectively. As usual we set 
j^(4) = ^ + j^j 4) with the standard choice [3, 4] 

J^j 4) = 8Gtr(0 t 0), J^j 4) = 2Htr[dld L + dld R ], (2) 

where fa = (qR) } a{q L ) l a, (<fc)oi = £abc£ijk(qL) J b C(qL) k c , and (d R ) ai = £ a bc£ijk{qR){C{qR) k c , 
with a, b, c and i, j, k the color and flavor indices, and C the charge conjugation operator. 
The flavor U(3) generators % a (a = 0, • ■ • , 8) are normalized so that tr[T a Tf,] = 28 a b, and 
Ta and X A i with A, A' = 2,5,7 are antisymmetric generators of flavor and SU(3) color, 
respectively, if ^ is invariant under SU(3)l x SU(3)^ x U(1)a x U(1)b symmetry. The 
interaction if^ 4 ^ produces attraction of qq pairs, leading to the formation of a chiral 

condensate. Similarly ifj 4 ^ leads to attraction of qq pairs in the color-anti-triplet and 
spin-parity ± channel, inducing a color- flavor locked (CFL) condensate [1]. We treat 
the two couplings G and H as independent parameters. 

The six-fermion interaction in our model consists of two parts, if ( 6 ) = if^ +if^ 
if^ is the standard Kobayashi-Maskawa-'t Hooft (KMT) interaction [18], 



^ = -8K(det(j) +h.c.) . (3) 

This interaction is not invariant under U(1)a symmetry, which accounts for the axial 
anomaly in QCD due to instantons. Consequently the mass of the r\' meson becomes 
larger than that of the other pseudoscalar octet Nambu-Goldstone (NG) bosons (n, rj , K) 
for positive value of K. On the other hand, the term (3) makes the chiral phase transition 
first-order as a function of T at /! = for the massless three-flavor limit [19]. 

As shown in [7], the instanton can induce a coupling between the chiral and diquark 
condensates through a new six-fermion term: 

^ 6 J=K'(tr[(dld L )(t>]+h.c.y (4) 

It is this term that is responsible for the aforementioned low temperature critical point. 
We assume K' > 0, so that qq pairs in the positive parity channel, (J^) = —(dp), are 
energetically favored. We keep K and K' as independent parameters [12]. 

The condensates favored by the interaction if ( 4 ) + if ( 6 ) are the flavor- symmetric 
chiral and diquark condensates in the spin-parity + channel, defined by 

(fa) = U/2)% (d Lai ) = -(d Rai ) = (s/2)8 ai . (5) 



TABLE 1. Two sets of parameters in the present three-flavor NJL model: 
The momentum cutoff is fixed at A = 602.3 MeV [4]. The dynamical quark 
mass M and the chiral condensate % at vacuum are also given. 





m q [MeV] 


GA 2 


HA 2 


KA 5 


M [MeV] 


X 1/3 [MeV] 




Set I 





1.926 


1.74 


12.36 


355.2 


-240.4 


Set II 


5.5 


1.918 


1.74 


12.36 


367.6 


-241.9 



Here the condensate order parameters are x and s. 

It is straightforward to derive the thermodynamic potential at the mean-field level [12] 



<A (2tt) 3 



~ 2T J 70l[81n(l + e -^)+ln(l + e -</ r ) 
where A is a momentum cutoff to regulate the vacuum energy, 

U{x,s) = 6Gx 2 + 3H\s\ 2 - AKx 3 - ^K'\s\ 2 x, 



(6) 



(7) 



is a constant term which is needed to cancel double counting of the interactions, and 

(4 = V / (v / ^ 2 + P 2 ±M) 2 + (2A) 2 , (of = \J{^M 2 + p 2 ±ii) 2 +& (8) 

are the dispersion relations for the quasi-quarks in the octet and singlet representations, 
with M and A the dynamical Dirac and Majorana masses, defined as 



1 



M = m q -A[G--Kx)X + -K'\s\\ A=-2[H--K' X \s\. 



1*. 



(9) 



These equations imply that ^ < is energetically favored for non-zero m q , while s is 
generally complex; the thermodynamic potential is a function of |s| 2 . 



PHASE STRUCTURE AND DISCUSSION 

The phase structures can be determined numerically by looking for the values of % and 
s that minimize the thermodynamic potential in Eq. (6). We follow the parameter choice 
of [4]. We show in Table 1 two sets of parameters we adopt, Set I and Set II respectively. 
We vary the strength of the chiral-diquark coupling (the K' term) by hand. We work in 
the flavor SU(3) limit, assuming m u = m d = m s = m q for simplicity. 

We show in Fig. 1 the phase structures for Set I (massless case) in the upper panel, 
and those for Set II (massive case) in the lower panel. Panels (a) and (b) show the results 
without and with the ^T'-term; in (b) we have taken K' = 4.2Kq with = 12.36/A 5 as 
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FIGURE 1. The phase structure in the (ju, r)-plane in the three-flavor NJL model without (a) and with 
(b) the K' term. The Upper and lower panels present the results in the massless case I and the massive 
case II respectively. Phase boundaries with a second-order transition are denoted by a single line and a 
first-order transition by a double line. The dashed-dot line at high T in case II shows the chiral crossover 
line, while the dotted line in (b) denotes the BEC-BCS crossover. See [12] for further details. 



a representative value. The phase diagrams contain a CFL phase with s ^ with U(l) 
baryon number broken, and other two phases both characterized with s = 0, a Nambu- 
Goldstone (NG) phase with % 7^ and a normal (NOR) phase with either % = (in case 
m q = 0) or x ~ (in case m q ^ 0). From the two figures in the panel (a), we see that 
the current quark mass leads to the critical point on the high temperature side of the 
first-order line of chiral phase transition [6]. The critical point moves downwards with 
increasing quark mass m q since it acts as an external symmetry breaking source of chiral 
symmetry breaking and thus smears the strength of the phase transition. 

The effect of nonvanishing K' can be seen by comparing (a) and (b). We indeed see 
that the low temperature critical point shows up at the other end of the line of the first 
order chiral phase transition. This is , as discussed in [7], because the ^'-term acts as an 
external field for %, which turns the first-order chiral phase transition into a crossover 
in the CFL phase where s ^ 0. Note that the CFL phase in the panel (b) accompanies a 

nonzero chiral condensate % ^ induced by the anomaly mixing term • 

The axial anomaly, for sufficiently large chiral-diquark coupling K', not only triggers 
the low T critical point, but also drives a BEC-BCS crossover in the CFL phase, as 
discussed in [17]. Within an NJL-type model such a BEC regime appears for sufficiently 
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FIGURE 2. The phase diagram in the (/i,/f')-plane at T = for massless quarks, with the NG and CFL 
phases. A dotted line separates the CFL phase into BCS- and BEC-like domains. The critical point and 
the critical end point are denoted by P and Q, respectively. 



large pairing attraction, H, in the ^-channel [14, 15, 16]. The novel feature here is 
that the axial anomaly helps to realize the BEC regime through its contribution to the 
effective qq coupling. This can be easily seen by extracting from Eq. (4) the dominant 
zero mode (fyj ~ SjjX/2) contribution to the quark-quark interaction 

~ lK'\x\tr[(d R -d L )\d R -d L )]-lK'\x\tr[(d R + d L y(d R + d L )}. (10) 

The first term increases the effective attraction between quarks in the + channel, 
while the second term is repulsive and suppresses the pairing. Thus when the chiral 
condensate is nonvanishing, as in the NG phase, the axial anomaly helps the formation 
of a diquark BEC condensate. 

In fact it is possible to show that at sufficiently large K' there are nine diquark bound 
states with mass Md(jU, T) < 2M(jU, T) where M(jU, T) is the dynamical Dirac mass of 
quarks at /i and T . Each diquark complex scalar has quark number ±2 so that it feels 
chemical potential 2/1. Thus when 2/1 hits Mj}(fi, T) from below a BEC condensate must 
start to form. Then the condition for the onset of a BEC approaching from the NG phase 
(NG-BEC boundary) is given by the condition [20, 14] 

2 J u=M D ( J u,7). (11) 

In order to see how the BEC domain in the CFL phase grows as a function of K', we 
show in Fig. 2 the phase diagram in the (n,K')-plane for massless quarks. The first- 
order line separating the CFL and NG phases for small K' eventually terminates at the 
critical point P. On the other hand, for K' sufficiently large, a BEC regime of bound 
diquarks appears across a second-order phase transition at a critical chemical potential 
/l = Md(h,0)/2 shown by solid line; the phase boundary meets the first-order line at 
the critical end point Q. A novel first-order transition from the BEC to BCS regimes 



appears between P and Q, with discontinuous changes of both the chiral and diquark 
condensates. 

In conclusion, the axial anomaly, by driving a coupling between the chiral and diquark 
condensates, plays an important role in the many body physics of QCD, making the 
phase diagram extremely rich. For one, we demonstrated that it can indeed produce the 
low temperature critical point between the hadronic phase and the color superconducting 
phase predicted by the previous GL analysis [7]. In addition, we have shown that the 
coupling helps the formation of a BEC of diquarks via increasing effective quark-quark 
attraction in the + channel. As a result, a BEC-BCS crossover or even the first order 
BEC-BCS transition can be realized in the CFL phase. 

Finally we note that very recently the extension of our analysis incorporating the 
effect of heavy strange quark mass was reported [21]. It still remains an important task 
to extend the analyses imposing the charge neutrality and /3 -equilibrium conditions. 

The numerical calculations were carried out on Altix3700 at YITP in Kyoto University. 
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